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Abstract 



We introduce a novel approach for the investigation of spin-wave excita- 
tions in itinerant ferromagnets. Our theory is based on a variational treat- 
ment of general multi-band Hubbard models which describe elements and 
compounds of transition metals. The magnon dispersion is derived approx- 
imately as the energy of a variational spin-wave state in the limit of large 
spatial dimensions. A numerical evaluation of our results is feasible for gen- 
eral multi-band models. As a first application we consider a model with two 
degenerate orbitals per lattice site. From our results we can conclude that 
spin- wave excitations in strong itinerant ferromagnets are very similar to those 
in ferromagnetic spin systems. 
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I. INTRODUCTION 



The most prominent examples for metals with a ferromagnetic order are the elements 
of the iron group, namely iron, cobalt and nickel. Although the magnetic behaviour of 
these materials is a well-known phenomenon, there are still many open questions in this 
field (for a general introduction, see e.g., Refs. ). It is generally accepted that the 
basic reason for ferromagnetic order is the interplay between the kinetic energy and the 
Coulomb-interaction of the electrons. Nevertheless, it is still a matter of debate which kind 
of minimal model must be used for the description of ferromagnetic materials. Whatever an 
appropriate Hamilton may be, from a theoretical point of view one would expect that it will 
be a hard analytical task to find even an approximate solution for such a real many-particle 
problem. 

The simplest theory, which gives an explanation for metallic ferromagnetism is the 
Hartree-Fock-Stoner theory!. A surprising result of this theory is the statement that fer- 
romagnetism occurs in any system provided that the product of Coulomb-interaction and 
the electrons' density of states exceeds a certain amount. This statement is the famous 
Stoner criterion. This criterion usually leads to surprisingly small critical values of the 
Coulomb-interaction for which ferromagnetism is predicted to occur. At first sight one may 
consider this posteriori justification of the Hartree-Fock theory, which certainly 

fails for stronger Coulomb interactions. However, in some simpler model systems like the 
one-band Hubbard model it is well known that ferromagnetism requires very large Coulomb- 
interactions if it exists at all. This is a definite contradiction to Hartree-Fock theory, which 
indicates that the whole Stoner picture may be inadequate for a thorough understanding of 
itinerant ferromagnetism. Similar objections could be raised against spin-density functional 
theory which, like Stoner-theory, is based on an effective one-particle description. Despite 
their conceptual shortcomings, these theories quite successfully describe some properties of 
the iron-group elements, i.e. the magnetic moment or the shapes of the multisheet Fermi 
surfaces^. Therefore, a competitive strong-coupling theory must meet these apparent suc- 
cesses of spin-density functional theory before it will be taken seriously. 

To this end, we recently introduced a variational treatment of multi-band Hubbard mod- 
els with a general class of Gutzwiller wave functionsifi These models allow the description 
of real materials, for example the elements of the iron group. As a first application we stud- 
ied the ferromagnetic transition in a two-band model. In contradiction to the Hartree-Fock 
theory we found that ferromagnetism requires quite large Coulomb interactions. In par- 
ticular, we demonstrated the decisive role of the intra-atomic exchange interaction, which 
is found to be irrelevant in the Hartree-Fock approach. The Gutzwiller theory shows that 
finite values of this exchange interaction are essential for metallic ferromagnetism. Based on 
these results, we suggested that the complex atomic Coulomb-interaction has to be taken 
seriously in theories on itinerant ferromagnetism. In particular, it appears to be essential to 
take into account exchange interactions which form local spins according to Hund's rule. We 
are presently calculating physical properties of the iron-group elements from our correlated 
electron approach!. First results for Nickel show that our method is able to resolve all major 
discrepancies between experiment and spin-density functional theory!. 

Experiments not only provide information about ground-state physics, but also yield 
insight into dynamical properties of materials. It is found that metallic and insulating fer- 
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romagnets behave similar with respect to their low-energy spin excitations. In both cases, 
inelastic neutron-scattering experiments show pronounced gapless spin-wave excitation£0. 
The understanding of these excitations is very important since they govern the magnetic 
phase transition at finite temperaturesEl. Theoretical methods which allow the determi- 
nation of spin-wave dispersions in Hubbard models are quite rare. It is obvious that a 
convincing description of spin excitations can only be obtained starting from a qualified 
theory for the ground-state. However, almost all earlier theories on spin-wave excitations 
in ferromagnets are based on effective one-particle theories. Some of these theories lead to 
surprisingly accurate results compared to experiments. For example, in Ref. |TT]] a random- 



phase approximation for iron and nickel was introduced. In other approaches, the spin-wave 
dispersion is determined via a mapping of the itinerant system to a ferromagnetic Heisenberg 
model (see, e.g. Ref. p9[]). 



In this work, we present a theory for spin-wave excitations in itinerant ferromagnets 
which is based on the variational ground states in Ref. 0. Our paper is organized as 
follows: In Sec. || we introduce the general class of multi-band Hubbard models and the 
corresponding Gutzwiller wave functions. Our general approach on the spin-wave problem 



is presented in Sec. |TI[ In Sec. |3 we evaluate the spin-wave dispersion for our general 
class of multi-band Hubbard models. Finally, we apply our results to a model with two 
degenerate bands in Sec. |V|. Short conclusions close our presentation. Technical details are 
deferred to four appendices. 



II. HAMILTONIAN AND VARIATIONAL WAVE FUNCTION 



A. Multi-band Hamiltonian 



In this paper we consider the following general class of multi-band Hubbard models, 

Here, c^^ creates an electron with combined spin-orbit index a = 1, . . . , 2N {N = 5 for 
3d electrons) at the lattice site i of a solid. For simplicity, we assume that the orbitals do 
not belong to the same representation of the respective point-symmetry group. For example, 
in cubic symmetry this means that there is only one set of s, p, Cg and t2g-orbitals. In this 
case, one-particle-states l^o), which respect the symmetry of the lattice, lead to vanishing 
non-diagonal local hopping-terms, i.e. 

(•^ole+^Vl^o) • (2) 

This relation simplifies the calculations in this paper, but there is no fundamental obstacle 
to extend our method to a more general case. 

We further assume that the atomic Hamiltonian 

H^■,.. = Y.'.K+ W-'-^-'-c+ cXc,3C^^ . (3) 

(T (T1,CT2,(T3,(T4 

is site-independent and readily diagonalized. 
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iJat = ^r^r , (4a) 
r 

mr = |r) (r| . (4b) 

Here, we introduced the eigenvalues Er and the eigenstates |r) of H^t- The diagonahzation 
of Hi,t is a standard exercise (see e.g., Ref. |T^). Knowledge of the states |r) means that 



we found their expansion 

\r) = J2Tj^r\I), (5) 

in the basis of the configuration states |/). In these states, a definite set of spin-orbit states 
a is occupied 

|/) = |o-i, as, ...) = c+ c+ ■ ■ ■ Ivacuum) (ai < as < ...). (6) 
For details about the notation, see Ref. 0, Sec. II. 

B. Gutzwiller-wave-function and diagrammatic evaluation 

In Ref. we proposed the following wave-function for a variational examination of the 
Hamiltonian (|I]): 

I^g) =Pg|$o) =nAG|$o) . (7) 

i 

Here, I'I'o) is any normalized single-particle product state and the local Gutzwiller projector 
Pi-Q is defined as 



i:G 



nAr = l + 5^(Ar-l)mr. (8) 



To simplify our notation, we suppress the spatial indices wherever a misunderstanding is im- 
possible, e.g., we omitted the index i on the rhs. of eq. (^. The real variational parameters 
Ar may have values between zero and one, where these two limits generate the ground-state 
both in the uncorrelated (i^at = 0) and the atomic limit [Hi = 0) of our Hamiltonian (|l]). 

In Ref. we showed that the expectation value of the Hamiltonian (|l|) can be evaluated 
for the wave-function (|^) in the limit of infinite spatial dimensions. In this section we only 
summarize the main ideas of the diagrammatic derivation which are important for our 
treatment of the spin-wave-problem in the next chapters. For all details we refer the reader 
to Ref. i. 

First, let us consider the norm of the wave-function (|^), 



(^G I ^g) = n ("^o "^o) . (9) 

i 

The square of the local Gutzwiller-projector can be written as 
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^■G = l+ E (10) 

/,7-'(|/|,|/'|>2) 

where we introduced the (local) Hartree-Fock-operators 

nH^ = n-r, (lla) 
nf = n,-n° (lib) 

for 1 = 1', and 

= [n n^^] niuh {J = In I'; I = J U h; I' = J U h) (11c) 

o-eJ 

= n n (lid) 

for / 7^ The operators c^^ (c^^) QlldD should be placed in an ascending (descending) 
order. An explicit expression for the coefficients Xi-ij' in ([T0|) is derived in Appendix (0). 

When we apply Wick's-Theorem to the right-hand-side of eq. (^, all terms are repre- 
sented by certain diagrams with lines 

i^.'^' = (<l>o|c+,c,.,,|<l>o> (12) 

and local vertices Xi-ijr. The special form of the operator P^q in ( |T0|) has two essential 
consequences for the structure of our diagrams. First, the definition of the Hartree-Fock 
operators together with eq. (H) guarantees that there are no local lines, i.e., 

^T' = • (13) 

Second, the constraint |/|, |/'| > 2 requires that at least four lines meet at every local vertex. 
When we evaluate the expectation values in 



Hi-atJ =^Er (rhi-r)^^ and (14a) 



we obtain diagrams, which contain one (for Hi-^t) or two (for Hi) external vertices for the 
lattice sites i (and j). If such a diagram possesses at least one internal vertex, we have 
lattice sites, which are connected by more than two lines. Such diagrams vanish in infinite 
dimensions and therefore we concluded in Ref. that the expectation values ([1^ only 
include diagrams without any internal vertex. Thus, we can write the expectation values 
{rhr)^^ = {rhi-r)^^ in (Pa]) as 



mr = (rhr)^ = Xlm'^ = (mr)o (15) 
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with uncorrelated expectation values (...)o = (•••)$„• This relation allows to replace the 
original variational parameters Ar by the new parameters mp. The expectation value for a 
hopping term in ( |14b| ) becomes 



r,r' 



(16a) 
(16b) 



X /i/i'v"'(^'u-)"'°'^rt(/ua)^(^'u-),r^r^,7'^/,r'. 
i,i'{(70,r) 



Here, the fermionic sign function 



/i=(/Ua|c+|/) 



(17) 



gives a minus (plus) sign if it takes an odd (even) number of anticommutations to shift the 
operator to its proper place in the sequence of electron creation operators in \lUa). Note 
that the numbers in ( |16a| ) are just the diagonal-elements of the matrix introduced in 



Ref. H, which is diagonal for our symmetry-restricted orbital basis (see Sec. |II A]) 



III. SPIN WAVES 



The theoretical examination of spin-wave excitations requires the analysis of the imag- 
inary part XT{q, E) of the transversal susceptibility0, which is given as the retarded two- 
particle Greenfunction 



GME) = ^{{SI,Sz))e 



L 



Here, we introduced the g-dependent spin-flip operators 

I 

= (Sty 



l,b 
-iqRl ; 



(18a) 
(18b) 

(19a) 
(19b) 



l,b 



in the Heisenberg-picture, where the sum includes all (L) lattice sites / and orbitals h. The 
magnetic excitations of the system are represented by poles of the Greenfunction GT{q,E) 
with energies E > 0. For our further analysis we expand the "spin- wave state" 



(20) 



in terms of exact energy-eigenstates 
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n 

The Lehmann-representation of (|18a]) , 



E-{E^-Eo)+i5' E + {En - Eo) + 15 



(21) 
(22) 

(23) 



shows that there are poles in Griq, E) for the energies En — Eq > with weights \Wn\ ■ 
In a ferromagnetic system the state | is also a ground state of H, since the operator 

S~_^ just flips a spin in the spin-multiplet of the ground-state l^&o)- Therefore, we can 

conclude that Gt{0,E) has one isolated pole for E — Eq = 0. Now we consider finite, 
but small values of q, and assume that the expansion (^T]) is still dominated by a narrow 
distribution of low-energy states. This scenario explains the pronounced peak in Xriq, E) 
for small values of E and |g |, which is seen in experiments and interpreted as a spin- wave 
excitation (see, e.g., Ref. [|10[)- Then, the spin- wave dispersion E^ can be identified as 



the position of this peak, and E^ is approximately determined by the first moment of the 

I 1 2 

distribution \W„\ , 



E, 



En En \Wn 
EnlWnf 



(^0 


H 




(^0 1 ^o) 



(^0 






(^0 







(^0 


H 




(^0 





(24a) 
(24b) 



It is still impossible to derive the spin-wave dispersion E^ from eq. (121) since we do not 
know the ground-state |\E'o) of our multi-band Hamiltonian (|l]). If we assume, however, that 
the variational wave function |\E'g) is a good approximation for the true ground-state |\E'o) 
we may substitute |\E'o) in eq. (|2^ by the variational wave function I^&g)- 
In the next section we will evaluate the "variational" spin-wave dispersion 



rpvar 
Eq 







^g) (^g 


H 


^g) 


(^G 


^^^^ 


^g) I ^g) 



(25) 



in the limit of large spatial dimensions. It should be noted that this quantity obviously obeys 
no strict upper-bound properties. Nevertheless, we expect that E^^ < E'"^'^ is fulfilled since 



the expectation values (p^) includes high-energy states which do not belong to the spin-wave 
excitation seen in experiments. 

In principle, transversal spin-excitations are given as peaks both in Xt(^; E) and 
Xriq-iEo — E) for energies E > Eq. In other word, we also had to consider the contri- 
butions from the Green function {{Sz^St))^ in our calculation. These contributions are 
identical to the second term in eq. ( p3D and we could include them by using the proper 
spin-wave state 
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(26) 



in our variational approach. However, the contributions from the second operator in (p6D 
vanish for q = and may be neglected for small values of |g |, where spin- wave excitations 
are observed in experiments. Nevertheless, there is no fundamental obstacle to extend our 
diagrammatic approach to the more general spin- wave state (P^. 



IV. VARIATIONAL SPIN- WAVE DISPERSION 



A. General considerations 



In order to determine the variational spin-wave dispersion (|25|) we need to examine the 
norm of the state \'^'^) = |\E'g) , 



(27) 



and the expectation values 

Hat 



^5 

9 



q 



Hi 



,iq{Ri-Rj) 



Strhk-rS:: 



q 



kl 



i,j,k,l 



(28a) 



:28b) 



Before we start to evaluate these quantities, it is necessary to discuss two general problems. 
First, let us consider the norm 



i,j,b,b' 



(29) 



in the special case q = 0, where S^_^ is just the total spin-flip operator S . When we assume 
that I^I/g) has the correct spin-symmetry, i.e. it is an eigenstate of 



S' = J2^^,^ and ^S'=(J2^'' 



(30) 



with eigenvalues Sq and Sq{Sq + 1), respectively, we obtain 

iV,-=o = 25g (^g I ^g) . 
Here, we used the well-known equation 

S+S- = S^- S'{S' - 1) 



(31) 



(32) 



for spin operators. In general, however, the wave functions {"^g), as defined in eq. (^, do 
not fulfill this symmetry. Therefore, it is necessary to introduce some additional constraints 
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^2 

on our variational parameters Ar in (H) to guarantee that \'^g) is an eigenstate of S" . In 
Appendix (0) we explain how these relations may be chosen. 

The second problem arises from the evaluation of g^dependent quantities in the limit 
of large spatial dimensions D. For example, let us consider the Hartree-Fock case, where 
I^g) = I'^'o) is a spin-polarized one-particle state with > We find 



1 ^ . 

X = ^ EE^"^'''"''^' (^0 |fit,TS..aSt.'aS>'.Tl ^o> (33a) 

6,6' i,j 

= E - <T<i + Z E e-(^-^^)i^2'^)'(^'T)p,(5-^)'(^'^A , (33b) 

b \ i^j J 

where, for simplicity, we assumed that the expectation-values P^J^ as defined in (|^) are 
diagonal with respect to the orbitals b, b' . For any finite value of q the sum in eq. ( p3b|) 



vanishes as 1/D. This means that the limits q ^ and D — > oo do not commute, because 
1™ ^ = E <T (1 - <i) ^ E «T - <i) = li^. -f- ■ (^'^^ 



To overcome this problem we will evaluate expressions like the sum in ( |33b|) using the realistic 



three dimensional band-structure of our Hamiltonian (0). This leads to results which are 
continuous in q and, consequently, reproduce the limit q = correctly . 

The norm of the state ( pT] ) for a correlated wave function |\1/g) 7^ \^o) will contain 
diagrams of an arbitrary order 1/D"'. In this paper we will only consider diagrams up to 
the leading order n = 1. At first sight one may wonder whether or not 1/D-terms need to 
be included in the ground-state energy expression as well. Fortunately, to order l/D, these 
diagrams only lead to a constant shift of the energy expectation-values for |\1/g) and |^^)- 
Thus, for our calculation we may neglect the 1 / D- contributions to the ground-state energy 
since we are only interested in the difference between these two energies. 

In the next subsection we evaluate the norm (|29|). The derivation of the expectation 
values (|28a|) and (|28b|) requires no additional techniques. The cumbersome calculations are 
deferred to appendices |^ and 0. 

B. Evaluation of the norm 

The norm of the state ( pTj) can be written as 

= E n ( (A;G^+5- A;g) P^,c)^ (35) 



?2 

$0 



When we use eq. (0) and apply Wick's theorem, we obtain diagrams with external vertices 
for the lattice sites i and j and internal vertices generated by the Hartree-Fock operators 
xi;ij'n]^fj, in PiQ. In Refs. P^,|5| it was shown that we only have to evaluate the connected 
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diagrams since the unconnected terms just give the norm of the Gutzwiller wave-function 
Nq = (\1/g I "^g)- For the evaluation of the first hne (|35|) we use the local relations 

5+5- = ^5+(r)mr , (36a) 
r 

5±(r) = Sir) [Sir) + 1] - s,{t) [s,{t) t 1] (36b) 

which follow from eq. (p2D. Here, we introduced the total spin 5'(r) and the spin-component 
Sz(T) of the atomic eigenstates |r). Hence, the expectation value of S^S~ is given as a linear 
function of the variational parameters and we may write ( pS] ) as 

" -^5_(r)mr (37a) 



LN, 



G 



+j E n { (A;g5.+A;g) {hcSiP,c) , (37b) 

where denotes the application of Wick's theorem and taking into account only the 

connected diagrams. 

For a further analysis of (|37b| ) we introduce indices D = (cTiCXi) for pairs of spin-orbit 



states CTj and the basic RPA-diagrams 

Pv'm = ptrS?) = -7 E e^'^^^-^^^pQ^^pjr ■ (38) 



L 



P^'{q) can be evaluated in momentum space as 



/-^ k k+q 

k 

with expectation values 

and a modified Kronecker-symbol 5^' = 6^^,. Note that, in contrast to the indices /, the 
order of the two spin-orbit states in 2) = (cridi) is significant. Here, its first and second 
element specify a particle which enters or leaves a vertex, respectively. 

For large spatial dimensions (i.e., up to the order 1/-D), the only contributions in ( P7b|) 
are RPA-type diagrams as shown in Fig ([^). The internal vertices generated by the 

operators P^q, which can be expressed in terms of Hartree-Fock operators xi-ij/fiffp (see 
eq. ([ToD ). For our RPA-diagrams we only have to consider two-fermion-operators in nffji- 
A general vertex with n incoming and outgoing lines is determined by the operator 

^ E E yj^hJUh {■■■ft'JUh,JUl2---Y<S,Q , 

/i,/2(/in/2=0) J{hMiJ) 
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with 

yjUh,JUl2 ~ ^JUll,JUl2JliJl2 ■ 

Thus, the internal vertex xS' = x^f^^'^^l at lattice site / stems from the term 

' -u (0-30-4) 

The round brackets indicate that Wick's theorem may not be applied to operators in different 
brackets. Thus, we can identify 

= ^(0-30-4) ^ f^t /oi ^(01,03), (02,04) • (43) 

An explicit expression for xiji is derived in Appendix |B[ 

Our infinite RPA-sum (see Fig ||) for the second term in ( |37b| ) leads to the following 
matrix-equation for fi§ (g*), 

n{q) = P{q) + P{q) ■ i ■ P{q) + P{q) ■ i ■ P{q) ■ 5: ■ P{q) + ... (44a) 
= P{q) + P{q) ■ X ■ n{q) , (44b) 

which has the solution 

f2(g) = (i-P(g)-x)-i-P(g). (45) 

Here, the dot "■ "indicates the usual product between two matrices. 

Finally, we have to examine the external vertices, generated by the spin-operators 
and S~ in ( P7b|) . First, we find for local operators 



PgS+Pg = J2 ^r'Ar |r') (r'| 5+ |r) (r| (46a) 
r',r 

= 5^Ar+Arv/:^|r+)(r| , (46b) 



r 

where the normalized spin-fiip-state 



was introduced. We may write eq. ( |46b|) as 



with 



PcS^Pg = J2 >^r,XrVS^)Yl Tj^^rT^jm^.j, (48) 
r h,i2 



rrij 



\{h, W (49b) 

oGJ agJ\(/U/') 
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and J = Ji n /2, /i = J U /, /2 = J U /'. 

The external vertices in our RPA- diagrams are met by just two lines. Therefore, we only 
need the following one-particle contribution of the operator rhi^ j^ in (^8]), 



(T1,(T2 

where we assumed that = I/2I. When we apply this result to the contribution of the 
external vertex ( ^7b| ), we only need to consider the first term with cxi 7^ cr2, because (Ti,(T2 
must have different spins. Thus, the expression for the external vertex, which stems from 
5",^ becomes 

^K..) = E Ar, Ar E ^/.r.^r^/.^tr^ (51a) 

r h,i2 

The expansion 

{...PoS-Po..]l^^ - Y: {■■■ i^tA.) •••};]* (52a) 

0-1,0-2 

= Y.{sU^.)*{-i^tAy-)l. ' (52b) 



01,0-2 



shows that the vertex-factor for the operator Sj in (|37b|) is given as 

s~ = (s^ V (53) 

(0-1CT2) \ {0201) J ' \ / 

When we consider and as components of vectors 5*+ and with respect to the 
indices T) we obtain the following final result for the norm 

= E ^- (r)mr + E St (g ) S^^ (54a) 
^ r J)i,D2 

= J2 S-{T)mr + S+- Q{q)- S~ . (54b) 



Note that the norm Nq in ( |54a|) will cancel out when we calculate the expectation values 

A similar derivation gives us the expectation values (^) of the kinetic energy and the 
Coulomb-interaction, see appendices and 0. There, the number of contributing diagrams 
is much larger (e.g., 18 for the kinetic energy) since we have up to four external vertices. 
Nevertheless, the numerical evaluation of these terms is a minor technical problem as soon as 
the wave-function I^E'g) has been determined by the minimization of our variational energy 
expression. First studies show that the application of our variational scheme to iron and 
nickel represents a solvable numerical taskQ. Further work in this direction is in progress. 
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V. APPLICATION TO A TWO-BAND MODEL 



In this chapter we will present the results for the spin-wave properties in a system with 
two degenerate e^-type orbitals per lattice site. The appearance of ferromagnetism in this 
model has already been discussed in Ref. 0, both for the Hartree-Fock and the Gutzwiller 
theory. Thus, we will only summarize these results here, before we start to consider the 
spin- wave dispersion. For all details, we refer to Ref. M. 



A. Ferromagnetic properties 

In a system with two degenerate orbitals the general atomic Hamiltonian (^) becomes 



(7,(J' 



cr 

In cubic symmetry the Coulomb- and exchange-integrals U, U', J and Jc are not independent 
from each other. Instead we have only two free parameters, because the relations J = Jq 
and U — U' = 2J hold. For the determination of the optimal variational wave function (^, 
we need the 16 eigenstates of the atomic Hamiltonian. The latter can be found in table I of 
Ref. i. 

In the one-particle Hamiltonian Hi we take into account first and second nearest neighbor 
hopping matrix elements. Furthermore, we apply the two-center approximation for the 
hopping matrix elements, which are chosen as t^^r ~ ^ ^'dlo- ~ 0-25 eV, and t'^^}^'^^ ■ t^dln"^^ ■ 
^d/^'' ~ l • '■ 0.1 (see Ref. [0 for the notation). The density of states for these 



parameters (see fig. 1 of Ref. 0) has a pronounced peak for the particle-density Tia ~ 0.3. 
For this band-filling we observe the strongest tendency to generate a ferromagnetic order, 
in qualitative agreement with the simple Stoner criterion. Therefore, we consider the spin- 
wave properties of our system only for this optimum band-filling. The numerical evaluation 
in Ref. did not respect the global spin-symmetry of the wave-function I^E'g)- In this 
work, we include the additional constraints on the variational-parameters Ar as described 
in Appendix 0, and we obtain almost the same state as in Ref. 0. In other words, in 
the two-band- model even the general variational ground-state I^^g); as defined in ([^, is a 
nearly perfect eigenstate of the global spin-operator. However, it is not clear so far, if this 
statement also holds for a more general multi-band model. 

In Ref. we found significant differences between the Hartree-Fock and the Gutzwiller 
theory. These differences have to be interpreted as a failure of the Hartree-Fock-theory since 
the variational space of the Hartree-Fock theory is included in our general class of Gutzwiller 
wave-functions. Although this statements holds for the corresponding RPA-theory as well, 
this theory is generally considered as the standard method in the context of spin-waves in 
itinerant ferromagnets.0 

The ferromagnetic phase diagrams for our model is shown in Fig. ^. In the Hartree-Fock 
theory ferromagnetism occurs for considerably smaller values of the correlation parame- 
ters. The most important difference lies in the role of the interatomic exchange J. In the 
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Gutzwiller theory a ferromagnetic ground-state exists only for finite values of J. This is 
completely different in the Hartree-Fock theory. There, the only relevant quantities are the 
Stoner-parameter / = ^^y^ and the density of states at the Fermi- level 'Do{Ep) which enter 
the Stoner criterion 

Wo{Ef) > 1 . (56) 

This means that even for J = finite values of U exist, where the system makes a transition 
into a ferromagnetic state. 

Another important difference between both theories occurs when we compare the con- 
densation energies -Econd, the differences between the energies in the paramagnetic and 
the ferromagnetic ground states. This quantity should have the same order of magnitude 
as the Curie-temperature Tq in itinerant ferromagnets. In Stoner theory we observe that 
-S'cond is typically of the order of U and therefore much larger than Tq. This is in agreement 
with the general observation that mean field methods overestimate the stability of magnetic 
order. On the other hand, in our correlated electron approach we find relatively small values 
for the condensation energy even for interaction parameters as large as twice the bandwidth. 



B. Spin- waves 

Fig. 1^ shows the spin- wave dispersion in (lOO)-direction for four different magnetizations. 
As the lattice constant of our simple-cubic lattice we chose a = 2.5A, the next neighbor 
distance in Nickel. Our Gutzwiller theory shows that the dispersion strongly depends on 
the magnetization, especially for small values of q. The line shows the respective fits for the 
function 

E,^=Dq\l-(3q') . (57) 

Note that experimental results of the region where quartic corrections become dominant 
are usually not very accurate, since Stoner excitations reduce the lifetime of spin-waves 
significantly. The values for the spin-wave stiffness D = 1.4eVA2 and D = 1.2eVA^ in both 
cases of almost fully polarized magnetizations, m = 0.26 and m = 0.28, respectively, are of 
the right order of magnitude for Nickel where D = OASeVA^. 

The spin-wave dispersion is almost isotropic, as can be seen for example in the inset of 
Fig. 1^, where the dispersion is shown in the directions (100) and (110). Such an isotropic 
behaviour was also observed in experiments on iron and nickel and it is actually somewhat 
surprising, since the band-structure in these materials is far from being isotropic. However, 
we should have in mind that even in a metallic system local moments are formed due to the 
electrons' correlations. Therefore, spin-excitations may be interpreted as spin-fluctuations 
in a system with localized spins. In such a system we have a generic isotropy when the 
exchange coupling is dominated by terms between nearest neighbors. 

Our results show that the magnetic excitations in strong itinerant ferromagnets behave 
very similar to those in systems of localized spins. These low-energy excitations are respon- 
sible for the magnetic phase transition which occurs for temperatures much smaller than 
the typical Fermi energies in itinerant electron systems. This observation is consistent with 
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the small condensation energy in our variational approach. To support this statement let us 
consider a Heisenberg-model 

Hs = — J ^ SiSj 

on a cubic lattice. In this system we have a spin- wave stiffness D = 2SJa'^. The value of the 
effective local moment in our itinerant system is given as 5* ~ 0.6 (see Ref. |^). Therefore 
we obtain J = ^ 0.16eV. For an estimate of the Curie-temperature we use the results 
from quantum Monte-Carlo calculations Tc = 1.44JS'^.li3 In this way we find Tc ~ 8 ■ lO^K 
which is the same order of magnitude as the condensation energy -Econd ~ 5- lO^i^. Thus, we 
can summarize that our variational approach gives a consistent picture of both of magnetic 
excitations and ground-state properties in strong itinerant ferromagnets. 

When the spin-waves are treated as non-interacting bosons, their contribution to the 
specific heat and the magnetization M(T) may be calculated from E^. The first order 
contribution stems from the quadratic term and is the well-known T^/^-law.lli The next order 
depends on the behaviour of for larger values of q. It is still not clear, wether or not 
the quartic term in (|57| ) describes the generic feature in the experiments. Some experiments 
indicate that the second term in M(T) is aT^, which could be explained with a linear 
behaviour in E^. Our dispersion E^ also allows to calculate such temperature-dependent 
quantities numerically. However, we did not analyze this for our two-band model, because 
these results could not be compared to experiments. 

VI. SUMMARY AND CONCLUSION 

In this paper we presented a variational method for the description of spin-wave excita- 
tions in itinerant ferromagnets. Our starting point was a general multi-band Hubbard-model 
which is an appropriate model for elements of the iron group. Earlier work showed that ferro- 
magnetism in metals requires substantial interaction strengths even in systems with orbital 
degeneracy. Therefore, we suspect that weak-coupling theories are not able to describe the 
physics of itinerant ferromagnets correctly. Our method is based on a variational study of 
multi-band Hubbard models with the help of generalized Gutzwiller wave functions. These 
wave functions yield the exact ground-state both in the uncorrelated and the atomic limit 
of the Hubbard model. Therefore, we expect them to describe reasonably the ground-state 
properties for finite values of the correlation parameters. 

From a theoretical point of view, spin waves are given as peaks in the imaginary part 
Xt{(1, E) of the transversal spin susceptibility . This quantity can be measured using inelastic 
neutron scattering. For q = the spin-symmetry of the ferromagnetic ground-state leads 
to an isolated peak 6{E) in xt- We assume that this peak is broadened only moderately 
for finite values of q. Then, the position of the peak, which is interpreted as the spin-wave 
dispersion E^, can be calculated as a static expectation-value of spin-operators in the ground 
state of the system. In our approximation we use the variational instead of the true ground 
state and calculate all expectation values in the limit of large spatial dimensions. 

Our results may be evaluated numerically for general multi-band models to describe iron 
and nickel. However, the main numerical problem is the determination of the optimum 
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variational wave-function for these realistic systems with a non-trivial atomic Hamiltonian. 
Work in this direction is still in progress. In this work we applied our method to a two-band 
model. Here, we found a behaviour which is in qualitative agreement with experimental re- 
sults for strong ferromagnets. The spin- wave dispersion Eg- is almost isotropic and quadratic 
for small values of q. For larger q we found quartic corrections, which are also seen in some 
of the experiments. The values of the spin-wave stiffness have the right order of magnitude 
compared to experiments. We concluded that the low-lying magnetic excitations in our 
correlated and itinerant electron system are similar to those in a localized spin systems. 
When we estimate the Curie-temperature Tc from our spin-wave properties and compare it 
to the condensation energy we find a consistent picture in our variational approach. The 
ferromagnetic phase transition is driven by spin-waves and the value for Tc is therefore much 
smaller than typical Fermi energies in itinerant electron systems. 
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APPENDIX A: GLOBAL SPIN-SYMMETRY 

The variational spin-wave dispersion (^) is gapless only if the variational wave-function 
|\1/g) is an eigenstate of the global spin-operator (see eq. (|30|)). This symmetry is not fulfilled 
for an arbitrary choice of the variational parameters Ar in (^. In this appendix, we present 
two ways to implement this symmetry. 

The first way to guarantee the global spin symmetry of I^&g) starts from the fact that the 
one-particle product-state |$o) in (0) i'^ general be chosen as an eigenstate both of S 
and S^. In this case, |^g) is an eigenstate of if the states |r) are chosen as eigenstates of 
the local spin-operator in z-direction. Then, the correct spin-symmetry of I^I/q) is ensured 
if |\1/g) is an eigenstate of S~^S~ . It is equivalent to demand that 

S^S-) =fi, (Al) 

since we may assume that the spin in I^&g) has a maximal component in z-direction. The 
left-hand-side of eq. ( [Al[ ) is just the norm for q = which was derived in ( |54a{ ). Thus, 



together with ( |54a| ) leads to one additional condition for the variational-parameters Ar. 



However, this condition is not very helpful, because equation (|54aD includes the optimum 



values Ap. This means that (|A1|) has to be included in the minimization algorithm, which 
determines the optimum wave-function |\I^g); numerically, this is a very difficult problem. 

In the following we propose a second, more feasible strategy for the implementation of 
the correct spin-symmetry. To this end, we arrange the orbitals on the atoms of our system 
into groups which carry the index of the respective representation D of the point-symmetry 
group. Then, we define the operators 

Nh=Yl "^-'ibs) , Md = J] mi-ibs) , (A2) 

i,b£D i,beD 
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for the gross and net number of electrons in orbitals of the representation D and with spin s. 
Again, we assume that each representation occurs only once. Under this condition, group 
theoretical arguments show that the following relation holds 

^mi.(bs) =^ni.ibs) - /£'(r)^';r ' (A3) 

6eD beD r(|r|>2) 

where 

/D(r) = E E l^r,7u(Mr • (A4) 

beD I[{b,s)ifl] 

Thus, we can write 

Mh = Nh- J2 fhi^Wr: (A5a) 

r(|r|>2) 

Mr = ^mi;r, (A5b) 

i 

and the Gutzwiller-projector (H) becomes 

D,s r 

Here, we already assumed that the parameters X(b.s) (= ^d) are the same for all orbitals, 
which belong to D. Further, we introduced 

Ar = Ar -11(^1.)"^^^'^ • (A7) 

D,s 

Now, we postulate the following conditions which ensure that |^I/g) is an eigenstate of the 
operator S : 

(I): Afj = Xjj, = As for all representations D,D'. 

(II): For all states l^g') , F^^V which belong to the same spin-multiplet (= F"^) with spin 



S, we have ApS^ = A = Ar^- 



For the prove of this statement, we can first conclude from (I) that the state |$o) is an 
eigenstate of 

D,s s 

where Ns is the number-operator for electrons with spin s. Thus, we have 

s r 

-S+S-JjAfl^o) , (A9b) 
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since 



for all s, r. We introduce the operator 



Ml 



(AlO) 



which has the property 



0. Then, condition (II) finally gives 



(All) 



such that |\1/g) is an eigenstate of S . 

Note, that condition (II) reduces the number of variational-parameters significantly, since 
all parameters Ar for states |r), which belong to the same spin-multiplet are now deter- 
mined by just one parameter Ar^. However, this restriction still allows different occupations 
of the several components, since 



also depends on the one-particle-state |$o). 
The one-particle-occupations 



D ; 



(A12) 



(A13) 



with b E D depends on the quantities nj^ = n(^bs) and m-^s^ (see eq. ([A3|)), 



m 



D 



n 



o- E /^(rfl^ri^ 

rs(|rs|>2)5,=-s 



(AM) 



However, the additional conditions (I) and (II) prevent us from the derivation of an analytical 
expression for mf,, since the parameters mps^ depend on m|) via ( [A7| ) and ( |A12| ). Therefore, 
the relation ( |A14| ) has to be implemented into our minimization algorithm with the help of 
appropriate Lagrange parameters. 



APPENDIX B: DIAGRAMS AND VERTICES 

1. The vertices and ^® ® ' 

In Section ( [TVD the vertices x®' (see eqs. (^)) have been derived in terms of the coeffi- 
cients xii>, which occur in the expansion (|lOl). Now we will derive an explicit expression for 
these coefficients. The operator ( [49 b| ) may be written as 
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rrij, 



h,l2 



n K + n (1 - -° - ^^') 

a€J'\{hUl2) 



= E 

J(/l,/2^J) 

Thus, Pq in (|1^) becomes 



JnJ' 



n< n (1--° 

crGj'Vj o-gJ'\{JU/lU/2} 



n 



HF 
J • 



^G-X^-^r ^J'u/i,r^r\j'u/2//i//i ( JI i _ ) 



J{h,l2^J) 



o-eJ'\J (TGJ'\{Ju/iUi'2} 
A comparison of the coefficients in ( P3|) and ([l0| ) gives 



n 



Ju/i,Ju/2 ■ 



xjuh,jui2- XI ^J'u/i,r^rtj'u/2.^/i./'/2 ( H i _ ) 



X 



JnJ' 



n-^ n 



creJ'vJ o-eJ'v{Ju/iU/2} 
In Appendix (0) we need an expression for vertices 

t'23'X>" _ P'{o-30-4)(o'50-6) 
^(o'i(J2) 

with three incoming and outgoing hues, which stem from a term hke 
Hence, these vertices are given as 



(T2U(T4 
0-6 



(o"3a"4)(o"5 0"6) /(T3 ^(TlUfTa /(T4 ^(T2U(T4 

(a"10"2) 



X(T3 ^CTiUfra ^0-4 !■ 
J G\ J (j^ J (J2 ^ <y(i 



((Tl ,0-3 ,0-5 ) (cr2 ,(T4 jCTs) 



(Bl) 
(B2) 



(B3) 



(B4) 



(B5) 
(B6) 
(B7) 



2. Diagrams 

In Appendix (0) and (0) we need some diagrams, which can be evaluated in momentum- 
space as follows: 



TI/('^10"2) — \ ^ iq(Rj~Ri) p(T10-4 po-3(T6 po-50-2 

'''^(<73(T4)((T5(T6)W ^ - 

^k+q k ^'T3^^3M^5<T4) ^a2^a2^(a3a4) 



fB8a) 



(B8b) 
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2. 



with 



0'30"4^(TiO-2 
k k+q 



-0"30-4 



E^ik{Ri-Rj)^a-iCT4, 



(B9) 



(BIO) 



3. 



(a'i(T2) 



(0-3 0-4) (0-50-6) 



iq{Ri—Rl)Mcrc,CFfi po-402 00103 
6 '^ij -f^ij 



E 



^0506 0402 0103 I ro-4 pK'^6)/'n ^ _i_ X03 pKo-6)/-.\ 

^fc +^02"'02^(oi03)lU j + V^oi^(a402)Wj 



(Ella) 
(Bllb) 



and 



{010-2) 



(03 04) (05 oe) 



= _ \ ^ Jq(Ri- Rj)f 0-5(76 pcr4a2 pa-ia3 

+ c:<2<S(^"') + ^:i^<i<:ol(o ) 



k+q 



(B12a) 
(B12b) 



7-r(oi(T2) l';?'! — \ ^ „jg(Rm--Ri)+'^l'^2 poi04 p0306 p0602 /'■R1Qq'\ 

^(a304)(0506)wJ - 2^ e r^^j- 2-i„, i-^i i^ij 



E 



0102^0104^0306^0502 _ ^06 V(<^504) 

k k k+q k '^3 03 (oi02)(oiO2) 



V'^oa Vl'^6)(0102)Wj '^o/o2'^oi'^02^(o306)WJ ''0i('03'^0i'^03^(0502)lU j + 

02 03 01 03 (0104)^ (0304) (0502)'^^'' (0302)'^^'' (0104) 



,0 ^0 771(0102), 



i^) - + (B13b) 



.0 ^0 77.(0102), 



APPENDIX C: EVALUATION OF THE ATOMIC INTERACTION 



According to eq. ( p^8a| ) we have to analyze the expectation values 



G 



^iqiR,-Rj) / ^ 



9 k ij(^k) 



S^'^k;vSj 



^g) +C.C. 



(Cla) 
(Clb) 
(Clc) 
(Cld) 
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This evaluation will be done separately for the three terms ( |Clb| )-( |CT^ ). 

In the first term (|Clb|) , we have to distinguish connected and unconnected diagrams. 
Here, the term "connected" means, that the lattice site k is connected to one of the lattice 
sites This condition necessarily requires that k is connected to both lattice sites i and j, 
because or generate a spin-flip. Such a process cannot be compensated in a diagram 
with only one of these operators, neither by the external vertex-operator rhk-r nor by one of 
the internal vertex-operators nffj/- 

The unconnected terms can be written as 



(C2) 



Pj^oSj Pj-G ] Pi-G 



When we ignore the restriction / 7^ fc, the sum over i and j gives just N^. Thus, we 
find the correct result for by substr acting all diagrams with i = k or j = k for an 



external vertex, or I = k for one of the internal vertices. These contributions are given as 



i= j = k 
i = k OT j = k 



I = k 



2NgS+ -niq)- s- 
NgS+ ■ n{q) ■ i ■ n{q) ■ s- 



Altogether, we obtain the following expression for the contribution of the unconnected dia- 
grams: 



{Clby 



Lmr - ^r^^ ( 25+ ■ ) ■ 5" 



mr' -|- ■ VL{q) ■ x ■ Cl{q) ■ 



(C3a) 
(C3b) 



where is given in (|54a|) . Note that only the second term m\{q) is relevant for our spin- 



wave-dispersion L^™^, since the first term Lmr is canceled by the respective ground-state 
contribution in eq. (^5|). 

For the connected diagrams in ( pib| ) we may distinguish between those diagrams with 
two or four lines, which enter or leave the external vertex. The vertex-factors with two lines 
can be evaluated from ([501), whereas the respective factor with four lines stems from the 
expansion 



m(o-io-2)(o-3o-4) 



ctA^) •••} 



*0 



,0-2 ,0-3,0-4 



E fai fa2 fas fa4''^'j(^ Y\. 



^"(0-1,0-2,03,04^/) 



og(oi,fT2,0-3,0-4) 



(C4) 
(C5) 



X 



[/o-i ^02 "Ju(oi,o-3)"JU 



"o-i"jUo-3"jUo-4 "o-3"jUo-i"jUo-2^ 



0-1 •'02 JU(oi,o-3) JU(o-2,04) 
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which is also derived for |Ji| = |/2|. This expression together with eq. 
following external vertex T§' for the operator 5*^^ 



r h,i2 



0-2) (0-3 0-4) 
2 



8|) leads to the 



(C6) 



whereas the corresponding factor for the operator S* is given as 



T 



(0-3 0-4) 

(t^lO-2) 



4^(0-40-3) 
((T2CTl) 



Using eqs. (|50| 



and 



/l,-f2 



(C7) 



we obtain the following expression for the vertices of the external-operator in (|Clb| ) with 
one or two incoming and outgoing lines, 



(0-10-2) 

,l2 



h,l2 



(C9a) 
(C9b) 



Il,l2 



— * 

When we define the vector M(T) of components M^(T), we may write the connected terms 
( pibp with z 7^ j as 



^2 



ml{q) 
ml{q) 



LN. 



G 



Si 



n{q) ■ s+ n{q) ■ s 



532 



2), 



LNg 



M(r) ■ f2(0 ) • i^T -niq) ■ S-+ T -niq) ■ 5^ 
S+ ■ Q{q) ■ M(r) ■ Q{q) ■ S' 



(ClOa) 
(ClOb) 

(ClOc) 
(ClOd) 



The tensors W!^^ ^^{q), f^^'^^ are defined in Appendix (Q. In fig. (||) all diagrams, which 
belong to the atomic interactions are presented. 

The connected diagrams in ( |Clb| ) with i = j are determined by the external vertex- 
operator 



(Cll) 
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This expression leads to the contribution 



^r(?"') = ^E^-(n MiV').^{Q).MiV) =m^(0 



For the evaluation of ( |Clc|) , we need the external vertex for the operator PqS' 
which is the same as the respective term in (^la|) for fixed F, 



A 



h,l2 



The only connected diagram in ( |(JlcD is therefore given as 



ml{q) 



LNq 



A+ ■ n{q) ■ S~ + c.c. 



Finally, we determine the contribution ( Cld|) as 



Here, we used the relation 

S~^mrS^ = S'+(F)mr+ • 
To summarize, the expectation value (|Cld|) for the atomic energy is given as 



H: 



at 



(C12) 
'rhrpG, 
(C13) 



(C14) 



(C15) 



(C16) 



(C17) 



APPENDIX D: EVALUATION OF THE ONE-PARTICLE ENERGY 

For the diagrammatic evaluation of the one-particle-energy we write the expectation 
value (128 b|) as 



Hi 



+ J2 ( e^^^(^-^') {SI 

i{^k,l) 



+ C.C 



c.c 



+ {Sk%a,SkCi..,)^^+C.C. 



(Dla) 

(Dlb) 

(Die) 
(Did) 
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Here, we already used the fact that the tight-binding parameters do not contain any lo- 
cal terms. In ( Dla|) we have to distinguish connected and unconnected diagrams. The 
unconnected contributions may be written as 



(Dla) 



:d2) 



X 



i^ii^k^l) m{jti,j,k,l) 
i m{y^i,k,l) 



Pj;GSj Pj-G ) P, 



>2 

m;G 



$0 



)2 

m;G 



Eq. ( |D2| ) can be evaluated, using the same arguments as discussed in connection with eq. 
^2f). This evaluation leads to 



where 



(D3a) 



(^S+ ■ Q{q) -S' + Y^ S-{T')mr' + S+ ■ ^iq) ■ x ■ Q{q) ■ S^" j (D3b) 



k=^l ak,<ri 

Before we start to discuss the connected diagrams, we should first consider the general 
structure of external vertices in (|Dla| ) at lattice sites k and I, where electrons are created or 
annihilated. For example, the vertex-function of the operator leads to 

{...Poc^Po...}l^ = J^ArAn J2 /.%':/WE^/i,rTr^,/4{-^/i-^.-}So (^5) 



r,r' 



Hence, our general problem is the calculation of vertex-contributions for operators rhj-^j^ 
with |Ji| — I/2I = 1. The first case is a vertex with only one incoming and no outgoing line. 



a' 



(D6a) 
(D6b) 

i'{u'ii') 

Note that eq. (|D5| ) together with (|16b|) yields directly the g-factor ( |16b|) , which occurs as 
the renormalization factor for hopping-processes. 

Furthermore we have the cases with two (three) incoming and one (two) outgoing lines. 
This processes lead to the vertices 



(D7a) 
(D7b) 



/'((T',cri,(72^/') 



<TG(cr',cri,CT2) 
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and 



(o"ia'2)((T30-4) / / 



/2 



cr 



ctG (cTl ,cr2 ,(J3 ,(T4 ,cr' J 



(D8a) 
(D8b) 



X 



-'^{o'l ,CT2 ,0-3 ,0-4 ,ct' 

icri io-2 "/U{cri,(73,cr')"/U(o-2,(T4), ^ cr' \^"cri "/U((T3 .cr') -fUcr4 " ai" I\j(as,u'f IU(T2 ) 



Using (p5|) we may now define the following vertex-functions for a single creation- 
operator as it occurs in ( pia| ) 

r,r' i{o-fi) h,i2 

r,r' /(o-^/) h,i2 
The respective vertices (5xi(c"; cr')? Qxid'(^! ^r annihilation-operators are given as 



Q 



(0-2 o-i) 



a, a 



Q{a^a2)(azai){,Cr,a') = Q 



(0-2 (Tl) (0-4 0-3) 



(a, a') 



(DlOa) 
(DlOb) 



Now we can determine the connected diagrams in ( pia ). First, we have the following 
terms for i j 



E vm2 e + 



(Dlla) 



X ^^i + x-n{q)j -S 



2)i 



i + x-n{q)) -s 



^'(^"') = ^ E vm2 E ^(o;o2)(o4.2)(o ) x.{i+m-i 



01,02 



X)i,D2,S3,S4 

©2 



AT, 



J34 



(Dllb) 



^'(^"') = ^ E E ^(o:o2)(oio2)(o ) 1+^-^(0) 



X)i, 1)2,2),, 1)4 



X)i,D2 



1 1)3 



©4 



(Dllc) 
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£%Q) = ^ Y1 5^ (<5i,(^l,^i)<5Dz(0^2,0r2) +(5i^((Ti,(T0<53.(o^2,C^2)) 



0"1,(T2 jCTjjO-j XlijXIz 



■ s- 



(Tl,(T2,0-j 



(Dlld) 



-I 



Dz 



(Dlle) 



1 „. _ „/ Si,Sz,S3,S,- ^ ^ ^ ' ^ ^' 



0"ljO'2,0"i 



X ^1 + ^](o) -x 

X [fi + x-fi(g)) -5 
AT. 



3, 



(Dllf) 



G 



^^)=iv: E V 

01,02,0-j^ 



■ E (^i. 



33, 



(Dllg) 



(ai, aOEj--) (0 ) + g^.x,, (ai, (0 ) 



-1(0201) 



Sz 



(Dllh) 



G 



AT.- 



X 



E v^E fe.(-i'-i)^(:;::)(o) +^^.(-1'-!)^^^^^^^^^^^^ 

01,02 ,oi Si ^ ^ 



(Dili) 



©1 



G 



E v/M^E^S 



01 02) (01 02) 



(0) 



3i 



X 



(Dllj) 



J 33i 



Expressions for the diagrams E^'^, U^^^^, V^^^^, and can be found in Appendix (0). 

Using eq. (|C11| ) we derive the connected diagrams in ( pia|) with i = j as 

e'^iq) = ^ E v/^EE^r'5-(nV^(^.2)(oio2)(0) + 5 -^(0)) -Mir')] , (D12a) 



e'%q) 



01,02 



Di r' 



i,o2,oi Si r' ^ ^ ^ ^ I 2 1; ^ 
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©1 



(D12b) 



For the evaluation of ( pib|) we need the vertex-function for the operators 

PgcJ-Pg = ^rAp v^Mn E /i^rt/r/u.,r'_ E ^/.r^r"^,/,^/./. , (D13a) 



r,r' 



/(<T^/) h,l2 



r,r' 



h,l2 



with one (or two) incoming or outgoing hues. These functions follow directly from eqs. 
( iDTal ) and (p8^ 



r(a,a') = E ArAr'^S^F) E /.'^^/Wl E ^/i.r^rv.D^^Alc^') , (D14a) 
r,r' ^(o-^-f) -fi.-fa 

r+(a, a') = E ArA^ V^S^F) E /i^^/u.^/.r E ^/.r^rv.f^/i,/.(^') ' (D14b) 
r,r' H'^ii) h,h 

R.u.A^,^') = E^rAr'A/54fO E /i^^/W^ E^/i,rTrt,,^&7V) , (D14c) 
r,r' H'T^) h,h 

K^<^2)i^' = E ^rX^'^/sJ^) E /X:/u.r/,r'_ E ^/.rTr^,/,^t/?V') • (D14d) 



r,r' 



Note that for the vertices of the remaining operators PqS^c^Pq and PqS^c^Pq the rules 
( piOp apply. Now we are able to summarize the contributions ( |Dlb| ) as 



.14 



0'l:0"2>0"i 

X \ (i + x-n{q)] -S'^^ +C.C. 



(D15a) 



.15 



cri,<T2,ij[ Tlx 



X p{q) ■ S-^^ +C.C , (D15b) 
^''(^"') = ^ E E K(^i,al)g.,.(a2,a^) +r(a2,a^)g+^(ai,ai)) i^fX^Ug) 

^ / f ^ V 1 2 y 



X 



) ■ J ^ + c.c. . 



(D15c) 



The contributions from eq. ( Die ) are 

Ng 



E (r+(ai,aO (r+(a2,a^))* + {r{a,,a[)r r{a,,a',)) E^^l^^J^{q) . (D16) 



o"l,cr2,(Tj,(T2 



Finally, we need the vertex-function for the operator PqS^c^S Pq in ( |Dld| ) with one out- 
going line 



27 



r,r' 



in (|Dld|) . This gives us the contribution from eq. ( pid| ), 



(Ti,(T2,crj 



Altogether we may write the expectation value for the kinetic energy ( |D1| ) as 



G 



Hi 



18 



c=l 
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FIGURES 



FIG. 1. Multi-orbital RPA-diagrams for the matrix il. 

FIG. 2. Phase diagram as a function of U and J for the Hartree-Fock-Stoner theory (HF) and 
the Gutzwiller wave function (GW) for (a) = 0.29 and (b) = 0.35; PM: paramagnet, FM: 
ferromagnet. 

FIG. 3. Variational spin wave dispersion in (lOO)-direction, for the two-band model; ng- = 0.29, 
J = 0.2U, and the values U/eV = 7.8, 10, 12, 13.6 correspond to m = 0.12, 0.20, 0.26, 0.28. The 
inset shows the Variational spin wave dispersion in (100) and (llO)-direction for m = 0.2 and 
m = 0.28 

FIG. 4. All contributing diagrams in the evaluation of the atomic energy. The diagrams 
and are defined diS^ = l — Qx,^ = l + xQ. Dotted and solid lines indicate the arguments 
and q, respectively. 



30 



